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Abstract: This paper describes and proves a canonical procedure to decouple perturba- 
tions and optimize their gauge around backgrounds with one non-homogeneous dimension, 
namely of co- homogeneity 1, while preserving locality in this dimension. Derivatively- 
gauged fields are shown to have a purely algebraic action; they can be decoupled from the 
other fields through gauge-invariant re-definitions; a potential for the other fields is gener- 
ated in the process; in the remaining action each gauge function eliminates one field without 
residual gauge. The procedure applies to spherically symmetric and to cosmological back- 
grounds in either General Relativity or gauge theories. The widely used "gauge invariant 
perturbation theory" is closely related. The supplied general proof elucidates the algebraic 
mechanism behind it as well as the method's domain of validity and its assumptions. 
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1. Introduction 

Consider perturbations around backgrounds with one non-homogeneous dimension (namely, 
with co- homogeneity one). This category encompasses two important and widely studied 
classes: spherically symmetric, static backgrounds where the radial coordinate r is the 
essential, non-homogeneous coordinate, and homogeneous cosmologies where time is the 
essential coordinate. 

Technically, the problem may be reduced through separation of variables, expansion 
in terms of harmonic functions over the homogeneous spaces and dimensional reduction, 
to a Id action invariant under certain gauge transformations. The issues are to optimally 
decouple the fields while choosing an appropriate gauge. 

Recently, a specific case which belongs to this class was solved, that of the negative 
mode around the Schwarzschild black hole background 0] . The analysis followed by post- 
poning any gauge-fixing; obtaining the action in terms of a "maximally general ansatz" 
(one which takes account of isometries yet is general enough to produce all equations from 
the variation of the action) ; using the power of the action formalism to transform the fields 
and decouple them into two sectors, algebraic and dynamic; and finally solving the gauge 
away, altogether reducing a problem of 5 fields and 2 gauge functions into a master equation 
for a single dynamic field. It was suggested there that the procedure could be generalized. 

In this paper we make good that suggestion of ||: we formulate the generalization 
in the form of a theorem and supply the proof, along with its conditions and limitations. 
We start in section || by setting up the problem and notation, and distilling the data for 
the general case. Then we choose a path to tread which leads us to the characterization 
of the algebraic fields in terms of derivatively gauged fields which we define, and to the 
formulation of theorem ffl, which is our central result. After stating the theorem we prove 
it. 

Our procedure overlaps with "gauge invariant perturbation theory" which is widely 
used in the literature. In section || we discuss its added value, and we provide further 
detail on the application to spherical symmetry and to cosmological backgrounds. 
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2. Theorem 



Set-up. Consider a Id quadratic action with a gauge symmetry. We denote the Id variable 
by x, and the fields and gauge fields respectively by the vector notation 

4> = 4> l {x) i = 1, . . . , np 

t = t a (x) a = l,...,n G (2.1) 

where np, n G are the number of (real) fields and gauge functions, respectively. 

We wish to generalize the procedure of [|]] to decouple the action. Since at this point 
it is not clear how to do that, we start by distilling the data for the general case. 

The most general action quadratic in fields which is at most quadratic in derivatives 

is 

S = S[<f>,<f>] = j dxC[4>,$\ 

£[<f), 4>] := 4>' T L 2 (x) ft + f Li ft + fL a <t> 

:= L 2>ij (x) ft'ft>' + L Mi (x) ft ft ' + L 0>ij (x) ft <jP (2.2) 

where S[4>, eft] is a condensed operator notation for a quadratic form which operates on 
(f) l (x) which the second line describes in vector notation in terms of (L2(x), L\{x), Lq(x)), 
a triplet of np x np x-dependent matrices and finally the third line spells out the vector 
notation. Lq{x) and L 2 {x) are symmetric, while since the action is defined up to the 
addition of d x (^(j) T Lcf^j where L is any symmetric matrix, L\(x) can be assumed anti- 
symmetric. 

The gauging is linear in the gauge functions, independent of fields and we assume it 
is at most linear in derivatives, which is indeed the case when analyzing perturbations in 
GR or gauge theory. It is given by 

50 = GZ:=G 1 Z f + G o ti (2.3) 

where the condensed operator notation of the first equality is exchanged by a vector nota- 
tion in terms of (Gi(x), Go(x)), a pair np x n G x-dependent matrices. 
Gauge invariance ties C and G by requiring 

= i<5£ = £[</>, G£] (2.4a) 

for all (j? , £ a , and hence in condensed notation, which will be expanded later 

CG = 0. (2.4b) 

Altogether the quadratic form S[4>, cf>] holds all the data, including that of the gauging 
which must be (a subset of) its symmetries. Our purpose is to decouple parts of the action 
as much as possible. In other words, we wish to "diagonalize" S, to bring it to a canonical 
form. Of course S (or C) is not an ordinary matrix whose entries are numbers, but rather 
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they are differential operators, or equivalently polynomials in d x . These observations 
lead us to an algebraic approach to the problem, namely to find the invariant properties of 
the quadratic form S. 

Examples such as Q show that the action contains a mutually algebraic sector, namely 
a subspace over which the restriction of C is algebraic and hence identical to the restriction 
of Lq. How could that subspace be characterized? The following definition is central 

Definition 1. The image of G\ in field space is called the derivatively-gauged (DG) fields. 

Notation. Log denotes the restriction of C to the DG- fields. 
We may now state the theorem 

Theorem 1. Given a Id quadratic action with a gauge symmetry, 2 the following canonical 
procedure can be applied to decouple it into two sectors, algebraic and dynamic, and 
completely eliminate the gauge: 

A. The derivatively gauged (DG) fields are mutually algebraic. 

B. The DG-fields decouple from the rest after a shift in their definition (assuming the 
algebraic Ldq is invertible), or more precisely a certain shift in the definition of 
coordinates on field space. 

C. The shifted field coordinates are gauge-invariant. 

D. The remaining action is algebraically gauged 3 and as such each gauge function can 
be used to eliminate a field ending up with a dynamic action in np — 2 no gauge- 
invariant fields (under the genericity assumption that both Gx, to be defined below, 
as well as G\ are non-degenerate). 



Proof. All of the algebraic information is contained in the explicit expansion of the gauge 
invariance identity ( |2.4| ) into powers of d x (integration by parts allows to transfer derivatives 
from (j) to £ and vice versa so only the total number of derivatives matters). It is instructive 
to assume the special case where the L matrices are x-independent. In this case one 
may Fourier transform the action and obtain C as a matrix whose entries are quadratic 
polynomials in k, the Fourier conjugate of x, and similarly the entries of G are linear 
polynomials in k. The gauge invariance identity fl2.4b|) 

can be expanded into powers of k 



from k° up to A; 3 yielding 4 /^-independent identities (between constant matrices!), which 
is all the information one has, and they turn out to exactly suffice for the proof. Here 



1 An algebraic structure consisting of "vectors" whose entries are elements of a ring, such as the polyno- 
mials, rather than a field (such as the real or complex numbers) is called a "module" . Accordingly, matrices 
with polynomial entries are linear transformations between modules. We are especially interested in the 
ring of polynomials in a single variable d x . This ring allows for the procedure of division with remainder 
and is therefore said to be an "Euclidean ring" . 

2 Where the action is at most quadratic in derivatives, and the gauging is at most linear. 

3 Namely, Gi = when restricted to this sector. 
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however we shall start directly with the most general case where the L matrices do depend 
on x. 

Item A. From the third (and highest) order of ( |2.4D in d x we find 

= ^6£ = <f/ T L 2 G 1 Z" + 0(d x 2 ) . (2.5a) 
Since the equation holds for all £ a we deduce that 

L 2 G 1 = (2.5b) 

with the interpretation that the DG-fields do not participate in L 2 , the kinetic part of the 
action. 

Moreover, since 5S = identically so does the second variation vanish 
+ 2 (Gi £") T L 2 (G + G[)e + (G 1 O t la <n f + 

+ o(d x 2 ) = 

= (Gtt'f L 1 G 1 t" + O(0 x 2 ) (2.6a) 



where in the last equality we used (2.5b). Since the equation holds for all £ a we deduce 
that 

GjL 1 G 1 = (2.6b) 

This means that once the action is restricted to the DG-fields not only are the 2-derivative 
terms absent, but also so are the first derivatives 

G 1 Z} = (G 1 ti) T L (G l (2.7) 

and hence, by definition, the DG-fields are mutually algebraic. 

Item B. We split the field space $ into the DG-subspace DG with coordinates 4>%q and an 
arbitrarily chosen complement X with coordinates <j? Xl namely 

<5> = DG@X (2.8) 

or equivalently 

4> i = m G , fx) ■ 

The action can be cast into the following form 

£[<t>X,4>x] = 4>DG L DG 4>DG + % 4>DG (£m4>x) + &X [4>X , 4>X ] (2.9) 

through a step which involves integration by parts and yields the definition of the mixed 
part to be 

C m = L 1 d x + (L + ^L[) . (2.10) 
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The types of indices of Log are Ldq = LoGap and we have Ldgo/3 = Lo a /3, f° r Cm they 
are C m = C m ar and finally for the remaining A-part of the action the types of indices are 

£-X = £-X rs- 

The square may be completed by defining the shifted field coordinates 



■>dg '■= <Pdg 



+ L D 1 G C m (/)x (2.11) 



where we used our assumption that Log is invertible to achieve a decoupling of the action 
into two sectors CI^dGi&dg] (algebraic) and £x[4>x, <j>x] (possibly dynamic) 

C[<p,4>] = Cdg^dgAdg] + £x[(f>x,<t>x] 
C-dg ■= Log 

£-x '■= £-x — £-m^ L D g C m (2.12) 

where C m denotes the adjoint operator of C m . 4 

Let us make two comments on the decoupling. First, the change of variables into the 
shifted fields ( |2.11 ) augmented by the identity on (ftx, namely (ftx '■= 4>x r is an invertible 



change of variables, or more precisely, it has unit Jacobian. Second, note that Cx receives 
an extra term in the process ACx = — L C m ( 2.12| ), 5 C m is at most linear in 



derivatives since the DG fields do not appear in the kinetic term (see the proof for item 
A above). Combining that with Ldq being algebraic we deduce that Cx still contains no 
more than 2 derivatives. 

One may wish for a more geometric interpretation of the decoupling, especially since the 
choice of the <px fields seemingly introduced some arbitrariness into the process. Geomet- 
rically decoupling is tantamount to rinding DG 1 - - a subspace which is the ^-complement 
to DG. Indeed DG 1 - is given by the equations <\>%q = (while DG is given by the equa- 
tions (/)■£ = 0). Thus the arbitrary choice of the complement X is replaced by the unique 
^-complement DG^. At first, our intuition with vector spaces may lead us to assume that 
any sub-space U has an S-comp lenient. However, in our case where the components of 
vectors and matrices are polynomials rather than numbers 1 we should exercise caution and 
indeed it is found that the existence of the complement depends on the invertibility of the 
restriction of S to U. 6 



Item C. Requiring gauge invariance of the action in its decoupled form ( 2.12j ) we find that 
the piece which includes §dg is 

= l - 8C = 4> T DG L D ^ 54> DG + ... (2.13) 



'By definition V<j>i, 4>2 f dx(C m cf>x) 4>i = f dx<f)i (£, 



t, 



>■>) 



5 Such a process of decoupling or eliminating algebraic fields is commonly called "integrating out", a 
name which originates from a path integral perspective. 

6 Indeed, while every subspace U still defines an S-orthogonal subspace U ± , the two do not necessary 
span the whole space. For example take the module H? (over the ring of integers) with the quadratic form 
S = [2 — 1; — 1 m] where m is any integer and take U to be the subspace generated by u = [1 0]. Now U ± 
is generated by u 1 - = [1 2]; However u and u do not span 'M? , and that can be traced back to the fact 
that Su = [2] is non-invertible. 
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Since the vanishing should hold for any (j)DG and since L D q is invertible we immediately 
obtain the required 

Hdg = • (2.14) 

Note that in this proof we used the gauge invariance equation (|2.4j ) at orders 0, 1 in d x , 
thus completing an all-order use of that identity. 

Item D. By definition cpx are algebraically gauged. If we denote by Gx = Gqx the gauging 
when restricted to DG^, we define the dynamic space D to be the gauge invariant space, 
namely the quotient 

D := DG ± /Im{G x ) (2.15) 

Cx is naturally defined on the quotient D since = \ 5Cx = £x[Gx £, 4>x], and we call 
it the dynamic action, to distinguish it from Cdg which is algebraic. 
Let us count the dimension of D 

dimD = dimX — dim Im(G'x) 

dimX = n F - dim Im(Gi) (2.16) 

Therefore, whenever Gi and Gx are non-degenerate we have dim Im(Gx) = dim Im(Gi) = 
uq and dimD = np — 2 no, as required. □ 

Relaxing the assumptions. During parts of the proof we made the following genericity 
assumptions: that Ljjq is invertible, and that Gi and Gx are non-degenerate. We are not 
interested here in trying to relax these assumptions, though it could very well be that not 
all of them are necessary. For instance, if G\ were degenerate it would mean that one of 
the gauge functions appeared only algebraically and as such it could serve to eliminate a 
field and the situation would reduce back to the case of non-degenerate G\. 

3. Discussion and applications 

Our central result is theorem |l] regarding Id quadratic actions, which include perturbations 
around backgrounds with essentially one dimension, such as static spherically symmetric 
backgrounds or homogeneous cosmological backgrounds, in systems with gauge - be it 
either General Relativity of gauge theory. The result contains several messages 

• The gauge can always be eliminated in favor of gauge-invariant variables, while keep- 
ing the action local (in x). 

• Each gauge function "shoots twice" - first, it is responsible for an algebraic field 
("constraint"), which decouples from the dynamic action, and second it is used to 
eliminate a field. 

• The algebraic fields are precisely the "derivatively-gauged" (DG) fields (see definition 
0). 



The proof follows step by step the expansion of the gauge variation identity (2.4) 
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Open questions. Several interesting directions for generalization remain open. Here we 
discussed the quadratic action, but it seems plausible that higher orders in perturbation 
theory can also be written in terms of the gauge-invariant fields. If we consider actions with 
more than 2 derivatives there seems to be no reason for a decoupled algebraic sector, but 
still there should be some canonical form, just like the one for a matrix over polynomials 
in a single variable. Finally, it is completely unclear at the moment how to proceed in 
dimension - to 2 essential dimensions and higher. 

"Gauge-invariant perturbation theory". Our procedure overlaps with the widely used 
"gauge-invariant perturbation theory" - the formulation of the equations solely in terms of 
gauge invariant fields. Moncrief (1974) found the advantages of gauge invariant variables 
for perturbations around spherically symmetric background using a variational principle 
(much like the action method in Q). The constraints were exhibited for general back- 
grounds and a detailed analysis was carried for the Schwarzschild background. Motivated 
by the analysis of collapse of a star Gerlach and Sengupta (1978) ||] applied it to pertur- 
bations around time-dependent homogeneous backgrounds and generalized by introducing 
the most general perturbed energy-momentum tensor. The method of Bardeen (1980) Q 
is mathematically equivalent is some respects but aims for cosmological applications. See 
H and references therein for a review of this method when applied to cosmological per- 
turbations, and see (6|, Q for more recent discussions of perturbations of the Schwarzschild 
black hole in the gauge-invariant formalism. 

The author does not pretend to be much of a historian, but let us attempt to discuss 
the value added by the current paper. Comments and essential relevant references are 
welcome. 

• We prove the procedure in full generality thereby clarifying the mechanism responsible 
for it, namely the algebraic properties of S, the quadratic form which is the action, 
and in particular the identities encompassed by ( p.4| ). 

In the literature it is commonly stated that the purpose of introducing the gauge 
invariant variables is to get rid of gauge ambiguity. 7 Indeed gauge invariant variables 
are gauge invariant. However, this argument is misleading since the non-trivial con- 
tent of the procedure is that there exists an invertible and local transformation (local 
in x) to these variables. To the best of the author's knowledge this fact is proven 
here for the first time in generality. 

• Our proof relies on the action being Id, or more precisely on demanding that the 
transformation into gauge invariant fields be local only in a single dimension. In 2d 

7 For example "the formulation removes the annoyance and drudgery of dealing with ambiguities due to 
infinitesimal coordinate transformations" U (section I) as well as "We shall remove the necessity of referring 
to such solutions [gauge modes - BK] at all by introducing a set of ....'gauge invariants' " (section IV). "To 
have genuine physical significance gauge independent quantities should be constructed from the variables 
naturally present in the problem, here the perturbations in the metric tensor and energy-momentum tensor" 
|^] (section 3). "In this review we develop a formalism ... [of] gauge invariant variables... [T]he physical 
interpretation of the results in unambiguous." (section 1). 
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(and higher) C, G become polynomials in 2 variables du, fi = 1, 2; there is no unique 
definition of DG- fields, and there is no general principle to guarantee an algebraic 
sector. This is an important demarkation of the domain of validity of the procedure, 
which again appears here for the first time to the best of the author's knowledge. 

• Our procedure recognizes the importance of the "derivatively-gauged fields" (see 
definition [l]). 

General properties of the action. So far we discussed the solution of Id actions. Before 
we turn to applications let us discuss some general properties of the construction of these 
actions for perturbations around co-homogeneity 1 backgrounds. 8 

A priori our procedure stresses the importance of using the action formalism in GR 
which in turn requires the use of a "maximally general ansatz" - an ansatz such that all 
the equations of motion can be derived from varying the action. However, with hindsight a 
"short-cut" ansatz is possible: it is enough to compute the action as a function of the fields 
appearing in an ansatz with all the uq DG fields as well as some np - 2 no subset of the 
other fields to represent the dynamic sector. This is allowed because by gauge invariance 
the last riQ algebraically-gauged fields do not appear in the action anyway, but we shall 
not use it in this paper. 

Perturbations around co-homogeneity 1 backgrounds in the variables x may be dimen- 
sionally reduced in a straight-forward way. The tensor indices of each field are decomposed 
into x and the rest. For instance, the metric perturbation tensor hn V decomposes into 
the x-tensor h xx , the x-vectors h x i and x-scalars h{j where i runs over the other coor- 
dinates. The dependence of the fields on the homogeneous coordinates is expanded in 
harmonic functions, substituted back into the action and the homogeneous coordinates are 
integrated over. Finally the gauge functions are also dimensionally reduced via reduction 
of indices and harmonic expansion. 

Continuous space-time isometries of the background become local (gauge) symmetries 
of the action while discrete isometries (which operate on the boundary) become global 
symmetries. In a gauge theory, on the other hand, all the symmetries of the background 
become global symmetries. Accordingly, the perturbations are charged under the isometries 
which they break. Since the quadratic action may include only interactions between fields 
in the same representation (more precisely, in a pair of complex conjugate representations), 
we recognize an important mechanism that decouples the action into sectors. 

We shall now apply these principles to two main cases of physical importance. 

3.1 Spherical symmetry; Cosmology 

Spherical symmetry. Here we shall describe the counting of the dynamic degrees of 
freedom for the 4d Schwarzschild black hole, in the terms of this paper, and postpone the 
detailed analysis and the rest of the results to a later publication (work in progress). 

Namely, backgrounds where all but one coordinate belong to a homogenous space such asR 3 orRxS 2 . 
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Sector 


"Even" 


"Odd" 


fields 

gauge functions 


h rr h r t, h>t% 
f& 4 1} tr(^) ftg> 

erH w 


/! (2) h (2) h i2) 
ri ti 'Hi 

e(2) 


np 
n G 


7 
3 


3 
1 


dim(L>) = np — 2uq 


1 


1 


DG, \DG\ = n G 


'Vr * fc r*t 


n ri 



Table 1: 4d Schwarzschild space-time: counting the perturbation fields as an application of the 
general procedure of section ||. The action decouples into "even" and "odd" sectors according to the 
charge under the parity isometry. We list the fields and gauge functions in each sector. Altogether 
we "predict" 1 dynamical even field (the field appearing in the Zerilli master equation) and 1 
dynamical odd field (the field appearing in the Regge- Wheeler master equation). Fields which are 
vectors with respect to the sphere (those which carry an i index) are decomposed into two families 
of spherical harmonics: one which is a gradient and another which is divergence-free, accordingly 
the fields are denoted by a (1) or (2) superscript. A similar notation holds for the tensor field hij. 
Finally we list the derivatively-gauged (and hence mutually algebraic) fields in each sector in the 
DG row. 

Working in standard r, t, Q, O = (6, <p) Schwarzschild coordinates the isometries are 
U(l)t x 50(3)q. The essential dimension of the background which we called x in the 
general case is r, and the homogeneous 3 space fibered over r is Rj x§^. 

We wish to dimensionally reduce the action to a Id action in r. The metric pertur- 
bation tensor is decomposed into h rr , h rt , h tt , h r i, h t i, hij where i,j run over 9, (p. A 
harmonic decomposition of the fields means a Fourier expansion in the t direction and 
spherical harmonic decomposition in the angular coordinates, thus transforming each field 
h(r,t,Q) -> J2u;,im h ^,im(. r )- The § au g e functions split as well ^ -> (f r , £ t , 

An important role is played by the ^2 parity isometry. Fields whose parity is (— )' 
are called "even" while those with (— )' +1 are called "odd". Accordingly the action decou- 
ples into an even part and an odd part. Table |l| summarizes the counting of fields and 
demonstrates the formula for dimD (item D of theorem H). 

Cosmology. In homogeneous cosmology the essential coordinate is t, the homogenous 
fiber is M 3 with its 3d Poincare isometry. The 3d momentum operator may be diagonalized 
on the perturbation space, yielding the wave-number 3-vector k. k breaks the 3d Poincare 
group down to the 2d Poincare group. The translation part of the latter is not broken, but 
the SO(2) rotations may be broken according to the SO(2) tensor type of the perturbation. 
Accordingly we distinguish three sectors which the action decouples into: scalar, vector and 
tensor. Note that while in the standard presentation the terms scalar, vector and tensor 
refer to 50(3) here they refer to the "little group" 50(2). Table ^ summarizes the counting 
of fields and demonstrates the formula for dim D (item D of theorem [l]) . 
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Sector 


scalar 


vector 


tensor 


fields 


hyp hif 2, •) h 2,z t r J 


h ■ h ■ 


hij 


gauge functions 




ti 




np 


4 


2*2 


1*2 


n G 


2 


1*2 




dim(.D) = np — 2n G 








2 


DG, \DG\ = n G 


htt h z t 


hu 





Table 2: Homogenous cosmology: counting perturbation fields as an application of the general 
procedure of section [| The action decouples into scalar, vector and tensor sectors according to 
the charge under the SO(2) "little group" isometry. We list the fields and gauge functions in each 
sector, remembering in 4d to associate 2 polarizations with each SO(2) vector or tensor. Altogether 
we "predict" 2 dynamical tensor fields (the 2 polarizations of the graviton). If the k vector is taken 
to lie in the z direction and i,j run over x, y. Finally we list the derivatively-gauged (and hence 
mutually algebraic) fields in each sector in the DG row. If we couple gravity to matter, then as 
long as the energy-momentum tensor of the matter is spatially isotropic, the matter is guaranteed 
to couple only to the 2 algebraic scalar modes. 
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